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To search for new evidence of the chiral p-wave superconducting domain in Sr2Ru04, we 
investigated the unconventional local transport characteristics of a microfabricated Sr2Ru04- 
Ru eutectic junction. We found that the anomalous hysteresis in voltage-current characteristics 
[as reported in H. Kambara et ai: Phys. Rev. Lett. 101 (2008) 267003.] is strongly affected by 
dc currents, but not by magnetic fields. This suggests that dc current acts as a driving force 
to move chiral p-wave domain walls; a domain wall trapped at a pinning potential is forced to 
shift to the next stable position, thereby forming a larger critical current path and resulting in 
the anomalous hysteresis. 

KEYWORDS: Sr2Ru04, p-wave superconducting junction, eutectic, critical current, chiral domain, hys- 
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1. Introduction 

Layered perovskite Sr2Ru04 (SRO)^^ is now widely 
believed to be one of rare examples of odd-parity spin- 
triplet pairing superconductors.^' Its superconducting 
order parameter is represented as the so-called chiral p- 
wave state with broken time-reversal symmetry, which 
is similar to the superfluid A-phase of ^He.*-* The pure 
SRO phase has a sharp superconducting transition tem- 
perature (Tc) of 1.5 K (1.5-K phase), while the SRO-Ru 
eutectic system shows a remarkable enhancement of Tc 
onset up to 3 K with its resistivity gradually decreasing 
to zero at 1.5 K. The superconductivity observed at 3 K 
is called the 3-K phase. The 3-K phase is regarded 
as the interface superconductivity in the SRO region be- 
tween SRO and Ru inclusions because its superconduct- 
ing volume fraction has been found to be considerably 
smaller than that of the pure 1.5-K phase from specific 
heat measurement.'''' Tunnel junction experiments using 
Ru interfaces®' showed a zero-bias conductance peak 
due to Andreev resonance, suggesting that the 3-K 
phase shows non-s-wave superconductivity. Recently, we 
have performed local transport measurements on micro- 
fabricated samples of the SRO-Ru eutectic system and 
observed the development of the superconducting link- 
age channels connecting Ru inclusions with decreasing 
temperature below 3 K.^^^ From the measurement of the 
critical current (/c) in the linkage channels, we found that 
the 3-K phase also has odd-parity pairing symmetry, sim- 
ilar to the 1.5-K phase. This shows that the SRO-Ru 
eutectic system includes naturally formed p-wave super- 
conducting junctions in itself depending on temperature; 
the 3-K phase region-normal-state region-3-K phase re- 
gion as the superconductor(S)-normal-metal(N)-S junc- 
tion at 3 K > T > 1.5 K, while the 3-K phase region-1.5- 



K phase region-3-K phase region as the S-S'-S junction 
at T < 1.5 K (Fig. 1). 



(a) T~3K (b)3K>T>1.5K (c) K1.5K 
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Fig. 1. (Color online) Schematic image of superconducting chan- 
nels in Sr2Ru04-Ru eutectic system: (a) nucleation of 3-K phase 
superconductivity around Ru inclusions at 3 K, (b) formation of 
weak links as S-N-S junctions at 3 K > T > 1.5 K, and (c) 
transformation of junctions to S-S'-S junctions at T < 1.5 K. 



We observed quite anomalous hysteresis in the 
voltage-current {V — I) and differential-resistance- 
current (dV/dl — I) characteristics of these junctions be- 
low 2 K.-'^^^ "Anomalous" hysteresis implies that it is dif- 
ferent from the usual hysteresis observed in usual Joseph- 
son junctions (JJs). That is, (i) voltage discontinuously 
decreases when the V — I curve switches to the next 
branch with increasing current, and therefore, (ii) the 
hysteresis loop shows the opposite direction compared to 
usual JJs.^^-* Thus, these features indicate the emergence 
of the internal degrees of freedom of the superconducting 
state. One possible explanation of the anomalous hys- 
teresis is due to the existence of chiral p-wave domains 
reflecting the chiral p-wave pairing of SRO. The exis- 
tence of chiral domains in SRO has been suggested by 
several experiments; the domain size has been estimated 
to be approximately 50-100 /im from Kerr effect mea- 
surement^^^ and approximately 1 um from a SRO-Cu-Pb 
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junction experiment.-'^'*) 

To investigate such anomalous transport properties in 
detail and examine chiral domains, we measured the lo- 
cal transport characteristics of microfabricated SRO-Ru 
eutectic samples, both in-plane and out-of-plane direc- 
tions. The advantage of milling a sample down to the 
micrometer scale by a focused ion beam (FIB) is that 
one can extract each conducting channel from a complex 
network*^) and detect anomalous features without aver- 
aging over the bulk sample. In this paper, we report that 
the anomalous hysteresis is strongly affected by dc cur- 
rents, but not affected by magnetic fields. This suggests a 
dc-current-driven variation, which is explained by the 
discontinuous motion of the chiral domain wall trapped 
at a pinning potential, which is discussed below. 

2. Experiment 

Eutectic crystals of SRO-Ru were grown in an infrared 
image furnace by the floating zone method. The sam- 
ple preparation is as follows. First, a thin platelet of a 
crystal was prepared by cleaving a crystal along the ah- 
plane and mechanically polishing it down to a thickness 
of a few tens of micrometers along the c-axis. Second, 
the platelet was roughly cut into a small piece of submil- 
limeter width in the a5-plane. Third, the small piece was 
glued on a SrTiOa substrate using epoxy adhesive. Last, 
four contacts for transport measurement were formed on 
the piece using an evaporated gold film. Then, the sam- 
ple was milled by an FIB to form a "microbridge" be- 
tween the voltage lead contacts (Fig. 2(a)). Since elec- 
tric fields concentrate in the microbridge (typical size 
is shown in Fig. 2(b)), the measured electric resistance 
is predominated by the electric resistance of the micro- 
bridge. Thus, we can study local transport properties by 
focusing on the microbridge from a complex network of 
linkage channels spread over the entire sample. To mea- 
sure transport characteristics along ///c, we formed a 
crank-shaped path by milling two alternate slits along 
the c-axis by an FIB (Figs. 2(c) and 2(e)) in a man- 
ner similar to the fabrication of intrinsic JJs of high-Tc 
cuprates.i''' -'^^^ We confirmed the existence of Ru inclu- 
sions in the microbridge region by FIB milling after the 
transport measurements were done (Fig. 2(f)). Figure 
2(f) shows that a piece of Ru inclusion was embedded 
slightly below the top surface at the center of the mi- 
crobridge. We observed about six pieces of Ru inclusions 
including small ones between the c-axis slits. The trans- 
port measurements were performed by dc and ac meth- 
ods. In the ac method, ac modulation with an amplitude 
of 30 ^Arms and a frequency of 977 Hz was used. The 
sample was cooled below 1 K using a ^He cryostat. Ex- 
ternal magnetic fields were applied parallel to the c-axis 
for evaluating the magnetic field effects. 

3. Results 

3.1 R — T curves of FIB-milled samples 

Figures 2(b) and 2(c) show schematic sample configu- 
rations for 1 1 jah (sample 1) and ///c (sample 2), respec- 
tively (note: sample 1 is the same as sample D used in ref. 
11). In sample 2, we note that the resistance of the ///c- 
part is high and that of the ///a6-part is low (^18%) 



at 4.2 K. In this study, we regard the configuration of 
sample 2 as ///c. Figure 2(d) shows the differential re- 
sistance at zero-bias current (K) versus temperature (T) 
curves of samples 1 and 2. If we assume that the overlap 
path length (0.7 /im) along the c-axis only contributes to 
the c-axis resistivity, the ratio of out-of-plane resistivity 
(pc) to in-plane resistivity (pab) of this material is esti- 
mated to be Pel Pah — 40 at 4.2 K. This is comparable to 
previously reported data.^-* 




a region milled from side 



Fig. 2. (Color online) (a) Schematic sample configuration for lo- 
cal transport measurements (not drawn to scale), and configu- 
rations of microbridges of (b) sample 1 for // /ab and (c) sample 
2 for I//c, respectively. Sample 2 was made from sample 1 by 
milling slits along the c-axis. (d) Zero-bias R vs. T for samples 1 
and 2. (e) Scanning ion microscope (SIM) image (10 X 10 Atm^) 
from the side view of sample 2. (f) SIM image from the top view 
obtained after the top and side surfaces were slightly milled. A 
piece of Ru inclusion (~ 3 ^tm length) was embedded at the 
center of the microbridge. 



3.2 Anomalous V ~ I and dV/dl — I characteristics 

Figures 3(a) and 3(b) show typical V — I and dV/dl—I 
characteristics of samples 1 (Figs. 4(a) and 4(b) in ref. 11) 
and 2 obtained by the dc and ac methods, respectively. 
The most peculiar feature of these systems is that the ab- 
solute values of voltage abruptly decrease at the thresh- 
olds /thi — ±3 mA and /th2 = ±4.3 mA for sample 1, 
and at /thi = ±3.1 mA and /th2 = ±4.7 mA for sample 2 
with increasing absolute values of dc currents, which are 
shown by bold arrows in the V — I characteristics of Figs. 
3(a) and 3(b). For the ac method, since the dV/dl — I 
peaks or steps correspond to the IcS of the superconduct- 
ing linkage channels, i*-* we can see that Ic becomes larger 
with decreasing current rather than increasing current. 
Furthermore, no negative dV/dl—I value means that the 
unusual voltage decrease is a switching phenomenon in 
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which no intermediate states exist. Thus, we found that 
the anomalous hysteresis is a common feature observed 
not only in the in-plane direction but also in the out-of- 
plane direction. The peaks in dV/dl — 1 oil j jc (sample 
2) are considerably sharper than those of I / /ab (sample 
1) reflecting a quasi-two-dimensional conductor; a junc- 
tion behaves like a tunnel along the c-axis; on the other 
hand, it behaves as a weak link in the a&-plane. We note 
that the reproducibility of the hysteresis after thermal 
cycling (warming up to 300 K and cooling again) was 
amazingly high up to around 10 cycles for both samples 
1 and 2. Furthermore, we note that the anomalous hys- 
tereses were observed for most of the samples of SRO-Ru. 
Thus, the anomalous hysteretic feature is not an inciden- 
tal phenomenon, but an intrinsic property. 




I (mA) I (mA) 

Fig. 3. (Color online) Typical V — / and dy/dZ — 7 characteristics 
with anomalous hysteresis of (a) sample 1 at T = 1.3 K and (b) 
sample 2 at T = 0.80 K. The dV/dl values were normalized to 
the zero-bias dV/dl at 4.2 K. The bold arrows in V — I char- 
acteristics denote the anomalous voltage changes. The open and 
filled symbols denote the different sweep directions, from zero to 
maximum (open black), maximum to minimum (solid red), and 
minimum to zero (open blue). 



3. 3 Dc current sweep width dependence of hysteresis 

Next, we show that the shape of the hysteresis strongly 
depends on the dc current sweep width. Figures 4(b) and 
4(d) show dV/dl — I characteristics for different sweep 
widths of samples 1 (A-D) and 2 (E-J), respectively. For 
sample 1, no hysteresis was observed in (A and B). On 
the other hand in (C), a narrow hysteresis, compared 
with that in (D), was observed. For sample 2, no hys- 
teresis was observed in (E) and a negative-current region 
in (F). On the other hand in (G), a narrow hysteresis, 
compared with those in (H and J), was observed. Note 
that a small difference between the hystereses of (G) and 
(H and J) appears when current is decreased. Also note 
that the shapes of the hystereses of (H and J) are al- 
most the same except for the dV/dl — I peak heights. 
The experimental results show that hysteresis width in- 
creases in a steplike manner when currents larger than 



/thS (shown by bold arrows in Figs. 4(a) and 4(c)) are 
applied. We note that the difference in the presence of 
hysteresis in (F), i.e., hysteresis in the positive-current 
region, whereas no hysteresis in the negative-current re- 
gion, is due to a small fluctuation in Ith- For the same 
reason, the curve (H) shows a large hysteresis, which is 
similar to (J) rather than to (G), though the current 
sweep width of (H) is slightly small compared with /th2 
in the V — I curve. 




-4 -2 2 4 
I (mA) 



Fig. 4. (Color online) Dc current sweep width dependence of hys- 
teresis for sample 1 ((a) and (b)) and sample 2 ((c) and (d)). (a), 
(c) V — I characteristics as a reference, (b) dV/dl — I charac- 
teristics of sample 1 for (A) ±2.6, (B) ±3.1, (C) ±3.7, and (D) 
±5.0 mA current sweeps, (d) dV/dl — I characteristics of sample 
2 for (E) ±2.7, (F) ±3.0, (G) ±4.0, (H) ±4.6, and (J) ±5.0 mA 
current sweeps. The vertical dashed lines represent the sweep 
ranges. The open and filled symbols denote the different sweep 
directions, which are the same as those in Fig. 3. The curves are 
offset by -0.5 units in (b) and by -6 units in (d) for clarity. 



3.4 Magnetic field effect on hysteresis 

To investigate the magnetic field effect, magnetic field 
cooling (EC) and zero-field cooling (ZFC) were carried 
out. If the anomalous hysteresis originates from a mag- 
netic vortex, the shape of the hysteresis curve would 
change between EC and ZFC. For EC, first, we applied a 
magnetic field larger than the lower critical field Hd (~ 
70 G^^)) parallel to the c-axis at 4.2 K (130 G for sam- 
ple 1 and 1000 G for sample 2), and cooled the samples 
below 1 K under the magnetic fields. Next, we switched 
off the fields and measured dV/dl — I at zero field. On 
the other hand, for ZFC, we cooled the samples under 
zero field and measured dV/dl — I at zero field. Figures 
5(a) and 5(b) (5(c) and 5(d)) show the dV/dl — I curves 
of FC and ZFC for sample 1 (2), respectively. There was 
no significant difference between the dV/dl — I curves 
of FC and ZFC, for both the samples. Figure 5(a) also 
shows the data obtained under magnetic fields applied 
after FC. In any case, when the applied magnetic fields 
are higher than Hd , magnetic vortices can penetrate the 



4 J. Phys. Soc. Jpn. 



Full Paper 



Author Name 



samples. If trapped vortices make Ic variation, a differ- 
ence in the shape of the hysteresis between the dVj dl — I 
curves of FC and ZFC should be expected. However, no 
difference is observed between the dV/dl — I curves ob- 
tained at zero field before and after the application of the 
magnetic fields. These results imply that the anomalous 
hysteresis is not caused by the magnetic vortex effect. In 
contrast to the dc current effect, magnetic field showed 
no significant effect on the anomalous hysteresis. 
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Fig. 5. (Color online) Normalized ciV/d/ — / curves obtained after 
(a) FC in 130 G parallel to c-axis and (b) ZFC for sample 1 
measured at 1.3 K, and those obtained after (c) FC in 1000 G 
parallel to c-axis and (d) ZFC for sample 2 measured at 0.5 K. In 
(a), the dV/dl — I curves obtained a.t H = 130, 350, and G after 
FC are shown together. The curves are offset by -0.5 units for 
clarity. The numbers denote the measured sequence. The open 
and filled symbols denote the different sweep directions, which 
are the same as those in Fig. 3. 



4. Discussion 

The experimental results are summarized as follows, 
(i) Anomalous hystereses appear for both current direc- 
tions (1 1 1 ah and 1/ /c). (ii) The width of the hystere- 
sis curves changes discontinuously depending on current 
sweep width, (iii) Magnetic fields do not affect the hys- 
teresis. 

Next, we analyze these experimental results by us- 
ing the chiral domain model proposed in our previous 
study. ^"'^^ Using this model, we assume that the angular 
momentum transfer from a given chiral state {px—ipy) to 
the opposite chiral state (px + ipy) by dc current induces 
domain wall (DW) motion. Thus, dc current causes an 
unusual Ic variation as a function of position because Ic 
is expected to be proportional to the cross-sectional area 
of spatially inhomogeneous linkage channels. If a DW 
exists, it is preferentially trapped at a defect at the be- 
ginning because a small overlapping area is energetically 
favorable for different order parameters of antiparallel 
domains. Thus, Ic would increase with increasing dc cur- 
rents and cause the anomalous hysteresis. We also note 



that there is no DW between parallel domains. Thus, 
the anomalous hysteresis would appear only between an- 
tiparallel domains. 

Regarding feature (ii), if we assume that /th orig- 
inates from a pinning potential due to, for example, 
a type of lattice defect, the well-reproducible /th af- 
ter thermal cycle is explainable. Assuming that apply- 
ing dc currents reduces effective pinning potential, the 
trapped DW can escape when sufficiently large dc cur- 
rents, which can overcome the pinning potential, are ap- 
plied. This is analogous to the depinning phenomenon 
of charge density waves with electric fields exceeding a 
certain threshold value. ^'^^ Using the simplest model in 
this case, we consider a tilted washboard potential as 
the pinning potential, although the pinning sites are not 
aligned periodically in an actual case. Then, the effec- 
tive potential UcS is given as a function of the position 
X as UcS ~ —Ucoskx — alx, where U is the amplitude 
of the potential corrugation and k and a arc numerical 
constants. 

We should consider the variation in the superconduct- 
ing condensation energy due to the variation in the cross- 
sectional area of the channel at the pinning sites. For 
antiparallel domains, since a small cross-sectional area is 
energetically favored, an increase in the cross-sectional 
area means a loss of the condensation energy. There- 
fore, the energy increase is proportional to the volume 



7r(r 



rg)^, assuming a simple "horn"-shaped channel 



(inset of Fig. 6(b)), where rg is the radius of the channel 
at the origin {x = 0), which has a minimum cross sec- 
tion, r is the radius of the channel away from the origin, 
and ^ is the coherence length of SRO. Assuming that r 
is proportional to x away from a; = in the simplest 
case, i.e., r = fix (/? is a numerical constant), UcS is 
rewritten as a function of a; (without a constant term) 
as UcS ~ -U cos kx — aIx+{nQ IIc{0)^ 13'^ x'^£_, where 
fioHdO)^ /2 is the condensation energy per unit volume. 
We have roughly estimated the order of magnitude of the 
3rd term, i.e., the condensation energy loss, to be ^0.1-1 
eV, where noHc{0) = 0.023 T;^) ,^ = 66 nm (in-plane), 
3.3 nm (out-of-plane) at T = K for the 1.5-K phase;"^) 
a; « 1 ^m; and /? ~ 1/20 from a rough estimation of 
the ratio of the cross sections. On the other hand, we 
estimated U to be ~0.1 eV for an order of magnitude 
assuming that U is comparable to the typical pinning 
potential of a magnetic vortex. Thus, these two terms 
are comparable. 

Figure 6(b) shows examples of Ucs variation under dc 
currents for certain values of parameters. If we assume 
viscous movement for the DW, it will be trapped at the 
local potential minima. Therefore, the DW moves to the 
next local potential minimum at /thi and /th2 with in- 
creasing dc currents from zero (process: O^-l— >-2). On the 
other hand, reducing dc currents from Jth2 to Jthi (pro- 
cess: 2— >-3) does not change the position of the DW. As 
a result, a larger Ic is observed during decreasing dc cur- 
rent process. An important feature is that the DW moves 
back to the origin around / due to the quadratic 
potential. Thus, the hysteresis loop starts again from a 
lower Ic even in the region of negative currents, which is 
similar to the case of positive currents, and repeats the 
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Fig. 6. (Color online) (a) Schematic V — I characteristics and (b) 
model of effective pinning potential of domain wall as a func- 
tion of position of channel under dc currents flowing conditions. 
The open and solid (red) circles denote the trapped positions 
at the local minima with increasing and decreasing currents, re- 
spectively. The dashed quadratic curve shows the condensation 
energy loss. Inset: model of chiral domain trapped in horn-shaped 
channel. IcO increases {I'^q, /"g) with the cross-sectional area of 
the channel. The hatched light blue parts of the cross section of 
the domain wall correspond to the volumes of the condensation 
energy loss. 



chiral domain is usually discussed only in the context of 
the in-plane direction theoretically. Further experimen- 
tal work is necessary to validate the existence of chiral 
domains or to determine whether other models can be 
applied to the SRO-Ru system. 

5. Conclusion 

In summary, we have observed anomalous trans- 
port characteristics in microfabricated samples of 
the Sr2Ru04-Ru eutectic system, both in-plane and 
out-of-plane directions. The anomalous hysteresis in 
differential-resistance and voltage-current characteristics 
are explained by the presence of chiral domains. The do- 
main wall, which is trapped at lattice defects, is driven by 
large currents beyond the pinning potential. The absence 
of the magnetic field effect on the hysteresis implies the 
coexistence of the Px ± ipy and Px{py) states in a micro- 
bridge sample with a size of the order of a few microme- 
ters. All the observed features suggest the emergence of 
the internal degrees of freedom of the chiral p-wave state. 
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Generally in a chiral domain, edge currents are ex- 
pected to flow along OWs.^^'^"^-* Thus, it is naturally ex- 
pected that magnetic fields would interact with the mag- 
netic moment of the edge currents. Feature (iii), however, 
seems to contradict this speculation. One possible reason 
for this inconsistency is that the actual chiral domain is 
not a simple Px ± ipy state, but the coexistence of the 
Px ± ipy and Px{py) states. Because of the translational 
symmetry breaking at sample edges, the Px state, which 
is expected for the 3-K phase symmetry,^'') seems to be 
stabilized near the sample edges. In fact, we observed a 
rough tendency for the ratio of the 3-K phase to the 1.5-K 
phase to increase with successive FIB processing (smaller 
than ^10 /im in sample size) in i? — T measurements for 
various samples. This suggests that the coexistence state 
should be taken into account as the microbridge becomes 
smaller. If the Px state forms along the edge, the effect 
of the chiral edge currents would be reduced by Meissner 
shielding current flowing along the outside of the px ± ipy 
state. As a result, the magnetic field effect on the anoma- 
lous hysteresis would weaken or disappear because of the 
mixing of the time-reversal symmetry conserving states. 
This speculation is consistent with experimental results 
in which no spontaneous supercurrents are observed us- 
ing a scanning Hall probe or a superconducting quantum 
interference device (SQUID). ^^'^^^ 

From feature (i) it is found that, in the chiral domain 
scenario, the DW can move not only in the in-plane direc- 
tion but also in the out-of-plane direction, although the 
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